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FORFWOJ  (U) 

(S)  This  document  contains  the  written  contributions  Raytheon  has  made 
during  1972  to  the  development  of  an  RF  Scattering  Model  for  the  Ivory  Coral 
phenomena.  This  information  is  presented  in  the  form  of  five  independent  re¬ 
search  papers  addressing  different  aspects  of  the  problem.  The  first  three  of 
these  papers  can  also  be  found  in  the  proceedings  of  the  Prairie  Smoke  RF- 
Sca tter i ng-Model  Workshop  and  the  remaining  two  will  be  published  in  the  pro¬ 
ceedings  of  the  Prairie  Smoke  III  Data  Review  Meeting^.  Such  proceedings  serve 
as  effective  vehicles  for  getting  research  results  to  the  modeling  segment  of 
the  community  as  socn  as  possible.  However,  nov?  that  the  year's  work  has  been 
completed,  these  closely  related  papers  can  best  be  read  together  and  for  this 
reason  are  being  published  here  as  a  set. 


G.  D.  Thome 


1  Proceedings  of  the  Prairie  Smoke  RF-Scattering-Model  Workshop,  held 
at  Stanford  Research  Institute,  18,  19  July  1972. 

2  Proceedings  of  the  Prairie  Smoke  III  Data  Review  Meeting,  held  at 
Stanford  Research  Institute,  31  November  -  1  December  1972. 
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ASPECT  SENSITIVITY:  SOME  CALCULATIONS  AND  OBSERVATIONS  (u) 

by 

G.  D.  Thome 

Raytheon  Company 
Sudbury,  Massachusetts 

ABSTRACT  (U) 

(S)  In  this  paper  a  refined  model  of  the  earth's  magnetic  field  is  used 
to  compute  where  in  space  the  RAM  radar  would  see  returns  from  the  heated 
volume  over  Platteville  if  the  scattered  were  strictly  field  aligned  and  highly 
aspect  sensitive.  These  calculations  are  then  compared  to  the  on-frequency 
target  locations  measured  by  RAM  on  three  days  and  close  agreement  is  found. 
Backscatter  at  the  transmitted  frequency  was  observed  only  from  those  positions 
in  space  where  the  radar  looked  strictly  normal  to  the  magnetic  field  and  this 
was  so  even  though  the  true  reflectioh  height  for  the  heater  frequency  on  these 
days  was  well  below  (10-30  kilometers)  the  height  where  the  orthogonality  con¬ 
ditions  were  met.  It  is  concluded  that  the  scatterers  responsible  for  on- 
frequency  RAM  backscatter  are  highly  aspect  sensitive  and  are  not  confined  to 
a  small  height  interval  about  the  heater  reflection  level.  This  contrasts  with 
Arecibo  results  at  essentially  the  same  radar  frequency  in  which  on-frequency 
backscatter  is  observed  well  off  of  perpendicularity  to  the  magnetic  field  (52") 
and  is  observed  over  only  a  narrow  height  interval  (less  than  a  kilometer)  about 
the  heater  reflection  level.  It  appears  that  RAM  and  Arecibo  are  not  seeing  the 

same  scatterers  and  this  should  be  borne  in  mind  when  the  absolute  scattering 
cross-sections  from  the  two  sites  are  compared. 
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I  INTRODUCTION  (U) 

(S)  It  has  been  established  that  electron  density  irregularities  can  be 
generated  by  illuminating  the  ionosphere  with  intense  HF  radio  waves  and  that 
these  irregularities  are  strong  enough  to  support  detectable  backscatter  at  HF, 
VHF  and  UHF  frequencies J This  paper  is  concerned  with  the  aspect  sensi¬ 
tivity  of  these  scatterers  at  UHF.  In  particular,  attention  will  be  confined 
to  the  aspect  sensitivity  of  the  on-frequency  backscatter  observed  by  the  RAM 
radar  at  435  MHz.  The  intent  is  to  show  that  on-frequency  backscatter  is 
observed  at  RAM  only  from  those  positions  within  the  neated  volume  where  the 
radar  looks  strictly  perpendicular  to  the  earth's  magnetic  field.  To  show  this 
the  beam  width  of  the  radar  must  be  small  compared  to  the  dimensions  of  the 
heated  volume  and  this  is  true  only  for  the  UHF  frequency  at  RAM.  As  of  this 
writing  only  on-frequency  UHF  data  is  available  from  RAM  and  consequently  no 
analysis  of  the  off-frequency  ("enhanced  plasma  line")  returns  is  given. 

II  MAGNETIC  FIELD  MODEL  (U) 

(S)  The  principal  components  of  this  analysis  are  a  magnetic  field  model, 
a  geometrical  calculation  showing  where  in  space  the  RAM  radar  looks  normal  to 
the  model  field,  and  experimental  data  showing  where  in  space  on-frequency  back¬ 
scatter  at  435  MHz  has  been  observed.  The  model  is  that  of  Hendricks  and  Cain4 
in  which  the  earth's  field  is  represented  by  a  99-term  spherical  harmonic  fit 
to  surface  magnetometer  measurements.  Since  the  number  of  terms  in  the  expan¬ 
sion  and  the  number  of  magnetometer  sites  are  finite  and  since  the  most  recent 
measurements  used  were  taken  in  1965,  the  model  is  not  perfect.  A  comparison 
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(S)  Cont'd. 

between  the  predicted  and  the  observed  field  has  been  made  for  the  Boulder 
Magnetic  Observatory  (about  40  km  from  Plattevil le) .  The  results  are  shown 
below  in  Table  I : 


TABLE  I.  COMPARISON  OF  PREDICTED  AND  OBSERVED  MAGNETIC 
FIELD  AT  THE  BOULDER  MAGNETIC  OBSERVATORY,  BOULDER,  COLORADO 


PREDICTED 

OBSERVED 

(by  Hendricks 
and  Cain) 

1971 

Average 

Typical  Variations 

Quiet  Day 

Storm 

Average 

Error 

Inclination  (Dip) 

68°  04' 

67°  33' 

3-6' 

10-20' 

+31  ' 

Decl i nation 

,  13°  34' 

13°  17' 

5-10' 

20-40' 

+17' 

Although  the  average  error  in  inclination  and  declination  is  less  than  a  degree, 
these  errors  are  significantly  greater  than  the  irregular  variations  typical  of 
quiet  days.  During  magnetic  storms  the  field  may  fluctuate  by  as  much  or  more 
than  the  average  model  error,  but  such  storms  are  infrequent  (of  order  once  a 
month).  The  accuracy  of  the  RAM  measurements  is  such  that  a  half-degree  error 
in  the  magnetic  field  model  would  produce  a  significant  discrepancy  between  pre¬ 
diction  and  observation.  Consequently  a  first-order  correction  has  been  made 
to  the  Hendricks  and  Cain  model  by  subtracting  31'  in  inclination  angle  and  17' 
in  declination  angle  everywhere  in  space.  This  forces  the  model  to  correctly 
predict  the  average  1971  magnetic  field  at  Foulder.  The  model  is  then  used  to 
predict  the  field  in  the  ionosphere  above  Platteville. 
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III  ASPECT  ANGLE  CALCULATIONS  (U) 

(S)  The  second  step  in  the  analysis  is  to  use  the  corrected  magnetic 
field  model  and  the  location  of  the  RAM  radar  to  compute  the  magnetic  aspect 
angle  for  points  in  the  vicinity  of  the  heated  volume.  The  radar  frequency 
and  the  elevation  angles  of  interest  are  high  enough  so  that  refraction  can 
be  neglected.  Under  these  conditions  the  calculation  of  magnetic  aspect 
angle  for  a  spherical  earth  can  be  conveniently  dene  using  a  geometrical  pro¬ 
cedure  developed  by  Millman5  and  his  approach  has  been  used  here.  A  fixed 
altitude  is  chosen,  defining  a  spherical  shell  around  the  earth.  The  magnetic 
aspect  angle  from  RAM  is  then  computed  at  closely  spaced  increments  of  latitude 
and  longitude  and  the  results  plotted  on  a  geographical  map  of  the  Platteville 
area.  Contours  are  drawn  showing  where  the  magnetic  aspect  angle  is  within 
0.1°  of  90°.  The  procedure  is  repeated  for  a  series  of  heights  and  in  this 
way  a  three-dimensional  picture  of  the  90°  magnetic  aspect  surface  in  space 
can  be  developed.  The  results  are  shown  in  Figure  1.  This  figure  shows  that 
ihe  surface  of  orthogonality  in  space  for  the  RAM  radar  is  a  dome  centered  at 
roughly  43°  north  latitude,  104°  west  longitude.  The  top  of  the  dome  is  at 
about  310  km  altitude.  Over  Platteville  the  surface  of  the  dome  is  at  about 
295  km  altitude  and  the  surface  slopes  to  the  south.  A  100  km  diameter  circle 
is  drawn  around  Platteville  to  show  about  how  big  the  heated  volume  is  in 
horizontal  dimension.  The  point  of  this  figure  is  that  if  the  scatterers  were 
highly  aspect  sensitive,  RAM  could  only  detect  that  fraction  of  the  scatterers 
lying  on  the  domed  surface  of  orthogonality,  no  matter  how  large  the  disturbed 
volume  actually  was.  If,  for  example,  scatterers  were  generated  at  all  alti¬ 
tudes  along  a  100  km  diameter  beam  above  Platteville,  RAM  would  only  see  those 
that  lay  in  a  thin  disk,  100  km  in  diameter  and  tilted  slightly  to  the  south, 

with  the  northernmost  edge  at  about  300  km  and  the  southernmost  edge  at  about 
290  km. 
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(S)  FIGURE  1  POSITIONS  IN  SPACE  AT  WHICH  THE  RAM  RADAR  CAN  LOOK 
PERPENDICULAR  TO  THE  EARTH'S  MAGNETIC  FIELD  (contours 
of  90°  ±  0.1°  magnetic  aspect  angle  as  a  function  of  altitude  viewed 
from  White  Sands)  (S) 
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IV  COMPARISONS  WITH  EXPERIMENTAL  DATA  (U) 

(S)  The  RAM  target  location  data  presented  by  Minkoff  has  been  plotted 
in  the  plane  containing  RAM,  Platteville,  and  the  center  of  the  earth.2  Figure 

2  shows  where  the  domed  surface  of  orthogonality  (labeled,  '90°  magnetic  aspect 
contour")  cuts  through  this  plane.  If  the  scatterers  were  highly  aspect  sen¬ 
sitive  and  if  the  magnetic  field  model  were  perfect,  RAM  targets  would  appear 
only  along  this  contour.  This  figure  has  been  drawn  to  scale  so  that  RAM  data 
points  (given  as  range,  elevation  angle  pairs)  can  be  plotted  directly.  Figure 

3  is  a  blow-up  of  the  region  above  Platteville  with  several  additional  pieces 
of  information  added.  First  of  all,  a  stippled  band  has  been  added  to  the  90° 
magnetic  aspect  contour  showing  how  far  this  contour  would  move  due  to  the  0.1° 
variation  in  dip  angle  typical  of  magnetically  quiet  days.  This  shows  that 
there  is  about  a  5  km  uncertainty  in  the  computed  position  of  the  90°  contour 
because  of  background  fluctuations  in  the  magnetic  field.  RAM  target  locations 
for  three  days  have  been  provided  by  the  Riverside  Research  Institute2,6  and 
are  plotted  on  this  figure  as  circles  joined  by  straight  line  segments.  Data 
for  October  12th,  13th  and  14th,  1971  is  labeled  12,  13  and  14  respectively. 

The  circles  for  each  day  show  tht  southern-most,  center  and  northern-most  por¬ 
tions  of  the  target.  The  accuracy  of  the  target  position  measurement  is  set 
primarily  by  the  accuracy  with  which  the  elevation  angle-  of  the  return  car  be 
read  from  range-time-intensity  displays.  I  estimate  that  this  can  be  done  to 

an  accuracy  of  about  a  tenth  of  the  radar  beam  width  and  I  have  made  the  radius 
of  the  circles  this  size.  The  third  type  of  data  shown  on  this  figure  is  the 
true  height  of  the  F-region  peak  at  the  time  the  RAM  observations  were  made.7 
These  heights  are  shown  as  horizontal  lines  labeled  12,  13  and  14  for  October 
12,  13  and  14,  1971 . 
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(S)  FIGURE  3.  EXPANDED  DRAWING  OF  THE  EXPERIMENTAL  GEOMETRY 

OVER  PLATTEVILLE  (U) 
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(S)  The  first  point  to  be  made  from  Figure  3  is  that  the  observed  returns 
lie  on  or  very  near  the  90°  contour  of  orthogonality  for  RAM  suggesting  that 
the  scatterers  are  in  fact  highly  aspect  sensitive.  By  "highly  aspect  sensi¬ 
tive",  I  mean  that  the  aspect  sensitivity  is  so  great  that  no  echoes  can  be 
detected  at  all  unless  the  radar  looks  strictly  normal  to  the  field.  The  second 
point  to  be  made  is  that  the  layer  peak  was  well  below  the  altitude  where  the 
RAM  targets  were  observed.  The  reflection  height  for  the  heater  was  presumably 
at  or  below  the  layer  peak  and  consequently  it  appears  that  the  scatterers  re¬ 
sponsible  for  the  RAM  echoes  are  well  above  (by  as  much  as  30  km)  the  reflec¬ 
tion  level  of  the  heater.  It  is  suggested  that  irregularities  generated  at 
the  heater  reflection  height  are  "mapped"  up  (and  presumably  down)  the  field 
for  at  least  a  few  tens  of  kilometers  where  they  can  be  observed  by  RAM. 

V  CONCLUSIONS  (U) 

(S)  Two  principal  conclusions  have  been  reached: 

(S)  1.  The  scatterers  responsible  for  the  RAM  on-frequency  returns  at 

435  MHz  are  highly  aspect  sensitive,  that  is,  the  scatterers 
must  be  viewed  at  90°  to  the  earth's  magnetic  field  to  be 
detected. 

(S)  2.  The  scatterers  are  not  restricted  to  a  thin  disk  at  the  heater 

reflection  height.  The  RAM  echoes  usually  come  from  the  top¬ 
side  ionosphere  whereas  the  heater  reflection  level  is  on  the 
bottomside.  It  appears  that  field-aligned  irregularities  are 
mapped  along  the  earth's  field  for  at  least  a  few  tens  of  kilo¬ 
meters  to  a  position  in  space  where  RAM  can  view  them  at  normal 
incidence . 
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AN  RF  SCATTERING  MODEL  FOR  A  RANDOM  DISTRIBUTION 
OF  IRREGULARITIES  (U) 

by 

B.  Rao  Fendyala  and  G.  D.  Thome 

Raytheon  Company 
Sudbury,  Massachusetts 

I  INTRODUCTION  (U) 

(U)  In  1956  a  theory  of  RF  scattering  from  field  aligned  irregularities 
was  published  by  H.  G.  Booker  and  applied  to  the  problem  of  modeling  the  radar 
scattering  from  ionospheric  irregularities  in  the  aurora.1  The  starting  point 
for  this  theory  was  an  ad  hoc  assumption  that  electron  density  irregularities 
exist  in  the  auroral  ionosphere,  that  these  irregularities  are  randomly  dis¬ 
tributed  in  space,  and  that  they  exhibit  an  anisotropic  Gaussian  autocorrela¬ 
tion  function  having  a  greater  scale  length  along  the  earth's  magnetic  field 
than  across  it.  No  attempt  was  made  to  ex  pi  air'  why  these  irregularities  exist 
or  why  they  are  Gaussianly  distributed.  What  wa ,  accomplished  was  to  work  out 
the  relationship  between  the  physical  properties  (size,  density,  anisotropy) 
of  such  irregularities  and  the  radar  properties.  This  made  it  possible  to 
use  radar  measurements  to  infer  the  physical  properties  of  auroral  irregularities, 
it  being  left  to  the  plasma  physicist  to  explain  why  irregularities  with  those 
properties  do  in  fact  exist. 

(S)  Similarly  in  this  paper,  we  will  use  the  Booker  theory  to  produce  a 
scattering  model  for  heater  generated  irregularities.  We  shall  find  that 
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satisfactory  f’t  to  the  exr*' rimental  data  can  be  had  assuming  that  tne  irre¬ 
gularities  are  highly  aspect  sensitive,  that  their  RMS  intensity  is  of  the  order 
of  l°o,  and  that  their  autocorrelation  function  transverse  to  the  field  exhibits 
two  dominant  scale  sizes,  0.1  meter  and  0.5  meter.  No  attempt  is  made  to 
reconcile  these  properties  with  the  plasma  physics  of  the  heater-ionosphere 
interaction  because  a  theory  for  the  generation  of  such  irregularities  is  not 
yet  available.  Until  it  is,  we  offer  the  scattering  model  described  below  as 
a  basis  for  evaluating  the  usefulness  of  potential  applications,  as  a  guide 
to  the  properties  the  plasma  physicists  must  explain,  and  as  a  guide  for 
planning  future  experiments. 


II  SCOPE  OF  THE  MODEL  (U) 

(S)  We  seek  a  scattering  model  which  can  be  used  to  predict  the  radar 
cross-section  over  bistatic  as  well  as  nonostatic  (backscatter)  paths.  However, 
olmost  all  of  the  quantitative  data  available  for  testing  the  model  has  been 
collected  over  backscatter  paths  and  consequently  we  shall  begin  (in  Section  3) 
by  deriving  a  relatively  simple  model  valid  for  backscatter  only.  We  shall 
also  assume  in  this  initial  model  that  the  autocorrelation  function  of  the 
irregularities,  transverse  to  the  magnetic  field,  is  a  simple  Gaussian.  When 
the  predictions  of  this  model  are  compared  (in  Section  4)  with  the  experimental 
data  we  will  find  that  a  significantly  better  fit  to  the  data  can  be  had  by 
adopting  an  autocorrelation  function  which  is  the  sum  of  two  Gaussians  of 
different  scale.  The  scattering  model  is  then  rederived  for  the  double  Gaussian 
autocorrelation  function  and  generalized  to  include  the  bistatic  situation.  The 
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result,  then,  is  a  scattering  model  which  is  believed  applicable  for  bistatic 
as  well  as  for  monostatic  geometries,  which  has  been  validated  by  comparison 
with  experimental  monostatic  data,  but  which  has  yet  to  be  validated  over 
bistatic  paths. 

(S)  Our  model  deals  only  with  the  on-frequency  scatter  wh'in  the  heater 
is  operating  at  or  below  fQf2  in  the  ordinary  mode  with  a  vertically  directed 
beam.  CW  heating  is  assumed  with  a  heater/antenna  combination  producing  the 
same  on-axis  density  in  the  ionosphere  as  does  the  Platteville  heater  at  full 
power. 


Ill  SCATTERING  MODEL  (U) 


(U)  The  scattering  model  presented  here  is  based  on  the  theory  of  weak 
scattering  (Born  approximation)  developed  by  Booker1  for  the  case  of  noniso¬ 
tropic  Gaussian  irregularities.  The  voluma  backscattering  coefficient  derived 
from  the  theory  is  given  by  Booker  as: 


a 


e 


2  (-8tt2T2)  (  -8rr2 

TLexPj—  | 


( L2-T2 ) s i 


1  n2i|»  | 


(1) 


The  symbols  in  equation  (1)  are  the  same  as  used  by  Booker.  It  is  necessary  to 
convert  the  volume  scattering  coefficient  o  to  a  radar  cross-section  oD  as 

measured  by  a  backscatter  system  in  order  to  compare  the  model  and  the  observa¬ 
tions. 
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(U)  Let  aR  be  the  cross  section  of  an  equivalent  conducting  sphere,  then 
the  peak  backscattered  signal  power  is 

p  _  PTcGo>2jR 

R  "  (4,)3R4  <2 

where  G0  is  the  effective  on-axis  gain  (assumed  to  be  the  same  for  transmitting 
and  receiving  antennas),  PTq  is  the  peak  power  of  the  transmitted  pulse  and  R 
is  the  range  to  the  scattering  target.  Now  consider  an  elemental  volume  dv  of 
the  scattering  target  with  a  scattering  coefficient  cy  Since  the  isotropic 

volume  scattering  coefficient  is  4rrae,  the  received  power  from  the  incremental 
volume  will  be 

P-rG^A^4TraDdv 

dP"  =  '  (4tt)3  R4  "  <3> 


Integrating  equation  (3)  over  the  scattering  volume  one  gets  for  the  total 
ceived  power 


R  =  (4,)2R4 


lPr 


G^cr^dv 


Equating  (2)  and  (4)  lends  to 


°R  =  4u 


/- 

J  Px. 


~2  °Bdv 
ToGo 


(U)  The  following  assumptions  have  been  made  of  the  scattering  volume  in  order 
to  evaluate  the  integral  in  equation  (5). 

(U)  1.  The  scattering  volume  containing  field  aligned  irregularities 

is  cylindrical ly  symmetric  about  the  magnetic  field,  and  is 
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centered  at  a  distance  D  below  the  point  on  the  field  where 
the  radar  achieves  perpendicularity  to  the  field  (Figure  1). 

(U)  2.  The  plasma  wavelength,  x^,  is  taken  as  constant  with  a  value 

corresponding  to  the  specular  point,  This  seems  valid  since 
the  height  interval  over  which  significant  scatter  comes  from 
is  small  compared  to  the  scale  height  of  the  electron  density 
in  the  F-renion. 

(S)  3.  (a N/M) ^  =  (a  N/N)2  exp 

2 

(a N/N)  variec.  with  x  and  y  because  the  heater  beam  is  of  finite 
width  and  varies  with  z  since  the  irregularities  weaken  as  one 
moves  above  or  below  the  heater  reflection  level.  Since  the 
interval  in  z  over  which  significant  scatter  comes  is  small  com¬ 
pared  to  H,  it  is  only  the  value  of  (a N/N)2  at  z=0  which  counts, 
so  we  can  write 

-  _  2  2 

(AN/N)2  =  (AN/N)2  exp  j-D2/H2J  expj-j^^—-  J 

where  (aN/N)q  is  the  peak  value  cc  (a  N/N)  which  occurs  at  the 
heater  reflection  level. 

2  2  2 

(U)  4.  sin  ’ji  =  L  /R  since  significant  scatter  obtains  only  for  small 

values  of  z/R,  i.e.,  near  specularity. 

(U)  5.  The  antenna  has  a  Gaussian  beam  of  half-width  B  radians  so  that 


G  =  G.  exp 

-  2  2 
zl+z 

o 

L  (BR)2  J 

(z+D )‘ 


2  +  2  +  2 
\ri  w^  r 
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(S)  FIGURE  1.  BACKSCATTER  GEOMETRY  FOR  THE  HEATER 

DIAGNOSTIC  RADAR  (S) 
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Thd  transmitted  pulse  has  a  Gaussian  half-width  P  so  that 
PT  =  PTo  exp  {  ~(x2/P2)} 

(U)  On  incorporating  the  above  assumptions,  the  volume  integral  in  equati 
(5)  becomes: 


v 
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dx  so  that 


(U)  Cont'd. 

The  integrals  are  all  of  the  form  J e 
8  4  "  00 

oR  =  -j’jv'27  (A N/N)q  T2L  exp  |  -D2/H2  j  exp  j  — —  j 

v^T 


8tt2T2} 


N 
\/T 


1  1 

1/2 

2  1 

1/2 

2  V(L2-T2) 

2  +  2 

•(BR)2  w 2 

(BR)2  xV 

- 

1/2 


Since  L>>T  for  high  aspect  sensitivity,  the  last  square  rcot  term  may  be  sim¬ 


plified  to 


(BR)' 


8n2  L2 

Tv 


1/2 


2v/2ttL 


XR 


1  + 


4  7T2  L2B2 


1/2 


The  second  term  in  the  root  is  much  less  than  1  which  means  that  it  is  the  as¬ 
pect  sensitivity  which  sets  the  vertical  thickness  of  the  effective  scattering 
volume  rather  than  the  vertical  beam  width  of  the  antenna.  Hence, 


4tt5W2RT2  X 


(AN/N): 


°R  " 


1+(W/P)2  1 /2  [l+2(W/BR)2  ] 


1/2 


exp]-P2/H2(exp 


-8tt2T2 


(6) 


(S)  Equation  (6)  can  be  used  to  predict  the  backscatter  cross  section  in  terms 
of 

P,  the  Gaussian  radius  of  the  pulse 

B,  the  Gaussian  radius  of  the  antenna  beam 

R,  the  range  from  radar  to  the  specular  point  in  the  scattering 
volume 

T,  the  Gaussian  autocorrelation  distance  for  density  fluctuations 
transverse  to  the  magnetic  field 
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X,  the  operating  wavelength 

D,  the  separation  along  the  field  between  the  specular  point  and 
the  heater  reflection  level 

_ H’  the  e-folding  distance  for  the  variation  of  (aN/N)2  along  z 

(AN/N)  ,  the  mean  square  electron  density  fluctuations  of  the  irre¬ 

gularities 

XN’  the  plasma  wavelength  corresponding  to  the  ambient  plasma  fre¬ 
quency  at  the  specular  point 

W,  the  Gaussian  radius  of  the  heated  volume  in  a  horizontal  plane 
at  the  reflection  level 

The  parameters  P  and  B  are  related  to  the  commonly  expressed  half¬ 
power  beam  width  and  '■jlse  length  as 


B  -  0.60  (HPBW)  and  P  =  0.564  (ct/2) 


(S)  The  measured  radar  cross  section  corresponds  to  a  scattering  volume 
which  in  general  represents  only  a  fraction  of  the  total  volume  of  the  heater 
generated  irregularities.  In  a  situation  where  (W/P)  and  (W/BR)  are  much  less 
than  unity  which  means  that  the  radar  pulse  encompasses  the  entire  volume  of 
the  irregularities,  the  measured  cross  section  stands  for  the  total  cross 
section  of  the  disturbance  and  it  is  given  as: 


5  2  2  (a  N/NP 
4it  krRT2A - j—2 

AN 


exp  J -D2/H2  f  exp 


-8tt2T2 


The  radar  cross  section  measurements  reported  by  Fialer,  Evans  and  Lomasney2 
using  SRI  SFCW  sounder  observations  are  in  the  form  of  The  VHF  and  UHF 
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croSj  sections  obtained  from  RAM  radar  observations  reported  by  Minkoff3,4,5 
which  represent  oR  have  been  converted  to  oT  to  standardize  with  the  SRI  mea- 
surements,  facilitating  an  easy  comparison  with  the  scattering  model. 

IV  COMPARISON  BETWEEN  MODEL  AND  OBSERVATIONS  (U) 

Iv-A  Single  Gaussian  Scale  Size  (U) : 

(S)  Using  the  equation  (7)  which  is  appropriate  to  a  single  Gaussian 
scale  size,  the  total  scattering  cross  section  is  computed  as  function  of 
frequency  in  the  range  10-435  MHz  for  different  values  of  T.  The  values 
adopted  for  the  various  parameters  are  as  follows: 

W  =  90  km 

R  =  900  km 

^  =  54  m 
D  =  0 

The  value  of  (a N/N) 2  is  adjusted  in  each  case  so  that  the  computed  cross  sec¬ 
tion  matches  the  observed  at  the  low  frequency  end.  Figure  2  presents  a 
comparison  between  the  observations  and  the  curves  based  on  the  scattering 
model.  The  curves  shown  are  parametric  in  T  in  the  range  0.1  to  1.0  m.  The 
required  values  for  (aN/N)rm$  fall  in  the  range  0.3  to  7%  with  higher  values 
going  with  lower  T.  The  curve  for  T  -  0.1  m  estimates  c.oss  sections  well 
above  the  observed  for  VHF  and  UHF  whereas  T  -  1 .0  m  provide  at  the  other 
extreme  cross  sections  far  below  the  observed  values.  The  curve  for  T  =  0.3  m 
seems  to  be  the  closest  fit  to  the  observations  of  all  the  four  curves  presented 
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(S)  FIGURE  2.  A  COMPARISON  BET"FT':  7’>!r  or  c,r  "VAT  IONS  AND  A  MODEL  BASED 

ON  GAUSSIAN  IRREGULARITIES  (S) 
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It  should,  however,  be  clear  frotn  the  comparison  that  the  shape  of  the  computed 
cur/e  is  such  that  no  single  scale  size  can  satisfactorily  reproduce  the 
observed  frequency  dependence  of  the  cross  section.  It  is,  therefore,  attempted, 
in  the  section  to  follow,  to  explain  the  observations  on  the  basis  of  a  two 
scale  size  model . 


I V . B  Sum  of  Two  Gaussian  Scale  Sizes  (U) : 

(U)  We  shall  treat  in  this  section  the  presence  of  two  Gaussian  scale 

sizes,  first  for  the  general  case  of  bistatic  configuration  and  then  reduce 

the  equations  to  the  backscatter  case  of  present  interest.  The  volume 

scatterirg  coefficient  for  a  bistatic  case  is  expressed  as: 

_  2  .  2 

2  v  Sin  x  (  ) 

°v  ~  |  A  c  /  e  |  ^  P  I  K  “j^i  ) »  K^-nii ) ,  (8) 


where  P  is  the  power  spectrum  of  (a  e/e) 

Let  us  assume  that  the  autocorrelation  function  of  A e/c  is  a  sum  of  two  Gaussian 
functions  with  different  scale  sizes,  i.e., 


P(x,y,z)  =  «  exo  -1  /2 


(x2+y2)  z 


- 

+(l-a)  exp 

"(x2+y2)  z2] 

-1/2 

2  +  ~2 

- 

L  t2  l2J 

(9) 


Fourier  transform  of  the  above  equation  gives 


p (M.f)  =  (2tt) 3/2  |  a  T2Li  exp |~  -1/2  |t2(£2+p2)+L2<;2  j 


+  (l-o)  T2L2  exp  -1/2  |  T 2(£2+r,2)+L2r2} 


(10) 
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Set  £  =  K(/2-^  ) ,  *»  =  K(m2-m1)  and  f  =  K(n2-  n1 ) 
and  substitute  |  A e/e | 2  =  ( x/an)4  (aN/N)2  to  obtain 


o  =  (2tt) 3//2  ,r  sin2) 


(AN/N)‘ 


j«T2L,  exp  [  -1/2 


+  K2T2  )2  +  K 


‘L1  <"2-"l)2}] 


+  (!-«)  T^L2 


expf-l/2  |k2T2(/2-^)2  +  K2T2  (nym,)2  +  K2L2  (r^-n,)2}  {(11) 

Let  us  define  a  coordinate  system  centered  at  the  scattering  volume  with  z  axis 
along  the  magnetic  field  and  y  axis  along  the  magnetic  east -as  shown  in  Figure 
3.  The  locations  of  the  transmitter  and  receiver  be  (x]  ,  y1  ,  z1 )  and  (x2,  y2>  z2) 
respectively.  Now  consider  a  unit  scattering  volume  in  the  disturbance  whose 
coordinates  are  (x,  y,  z).  Then  we  have  for  the  direction  cosines 


L  .  1 

Ri 


(y  -  y} ) 


(z  -  Z-j) 


t  .  (x2  -  x) 
2  R2 


(y2  -  y) 


(z9  -  z) 


where  R-|  and  R2  are  the  ranges  from  the  transmitter  and  the  receiver  to  the 
scattering  volume. 
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(U)  FIGURE  3.  SCATTERING  GEOMETRY  FOR  BISTATIC 

CONFIGURATION  (U) 
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Let  us  denote 


(m2-m1 )' 

t 

(n2-n1 )' 


X-i  Xo  1  1 

+  -I  .  X(-  +  -  ) 


R1  R2 


R2  R1 


=  [A,  +  BxJ 


■  LA2  +  B*] 

'  [S3  +Bz] 


where  A-j ,  A2,  A^  -’’d  B  are  all  constants. 

Substitute  (12)  in  equation  (11)  and  integrate  over  all  space  to  get 


oc  Oo  00 


aRB  =  4 


■  W\l 

-  Oo  .CO  -  JC  ll- 


^  (j^  ■j  »mi  *  Hi )  G(  -/2  ,  -m2  *  —n2 ) 


avdx  dy  dz 


In  the  above  equation  the  transmitting  and  receiving  antennas  are  assumed  to 
have  identical  gain  patterns.  The  magnetic  aspect  in  a  bistatic  configuration 
is  given  from  the  direction  of  the  bisector  of  the  incident  and  scatter  direc¬ 
tions.  The  angle  <p  which  is  the  complement  of  the  angle  between  the  bisector 
and  the  magnetic  field  is  yiven  by  the  relation 

r  1/2 

sin*  =  (n2-n1)/[(i2-^1)2  +  (m2-m1)2  +  (n2-ni)2  ]  (1 

For  the  backscatter  case  of  present  interest,  we  substitute 

(^-|  »m1  ,n1 )  =  -(jp2,m2,n2)  =  (£,m,n) 

n  =  sin  ip  and  (/2+m2)  =  cos2i^  to  get 
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6 


=  (2tt)3/2tt2(aN/N)2  jaT^  exp  [—^1]  exp  (Lf -if ) sln^ J 


-8t 


+  0-«)  T2L2  exp|^-^ijexp^~(L2-T2)  sin^J  J 


(15) 


(S)  Following  the  steps  described  in  Section  3,  the  total  cross  section  is  given 
as : 


4tt5w2Rx  (aN/N)2 

aj=  - a 

an 


L  exp  j  -D2/H2||  aT2  eyp  y-1  J  +  (1-«)t|  exp[  “7^]  j (16) 


-8tt2T2 


From  equation  (16)  it  should  be  obvious  that  any  number  of  scale  sizes  can  be 
added  without  altering  the  form  of  the  equation.  Judging  from  the  results  pre¬ 
sented  in  Figure  2 ,  it  seems  that  a  combination  of  the  scale  sizes  with  T=0.1 
and  0.5  m  is  best  suited  to  obtain  a  reasonable  tit  to  the  experimental  measure¬ 
ments.  Figure  4  shows  a  fit  of  the  two  scale  size  model  with  the  measurements, 
and  in  addition  curves  are  presented  there  also  for  individual  scale  sizes 
for  comparison.  A  value  of  about  0.8%  for  (aN/N)rms  with  the  larger  irregulari¬ 
ties  nearly  twice  as  strong  as  the  smaller  was  found  necessary  to  obtain  a 
match  between  the  model  and  the  observations.  The  comparison  shown  here  makes 
it  clear  that  it  is  essential  to  allow  more  than  single  scale  size  to  exist  in 
the  scattering  volume  in  order  that  the  model  predicts  the  same  shape  as  the 
observations  seem  to  indicate  for  the  frequency  dependence  of  the  scattering 
cross  section. 
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(S)  FIGURE  4.  A  FIT  OF  THE  TWO-SCALE-SIZE  SCATTERING  MODEL  TO 

THE  OBSERVATIONS  (S) 
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RF  SCATTERING  MODEL  SUMMARY:  JULY  1972  (U) 

by 

G.  D.  Thome 

Raytheon  Company 
Sudbury,  Massachusetts 

ABSTRACT 

(S)  On  July  18-19,  1972  an  RF  Scattering  Model  Workshop  was  held  at  the 
Stanford  Research  Institute  in  California.  The  purpose  of  this  workshop  was 
to  arrive  at  a  model  of  the  scattering  phenomenon  which  is  consistent  with 
existing  experimental  data,  which  can  be  used  to  predict  the  performance  of 
hypothetical  scatter  communication  systems,  and  which  will  serve  as  a  focal 
point  for  planning  future  theoretical  and  experimental  efforts  in  this  area. 

We  have  produced  such  a  model  and  we  present  it  here  for  use  by  the  community. 
This  model  should  be  viewed  as  a  statement  of  the  Dresent  state-of-the-art 
rather  than  a  finished  product. 
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I  INTRODUCTION  (U) 

(S)  The  model  is  intended  primarily  for  use  in  predicting  the  performance 
of  communication  systems  using  the  Ivory  Coral  scattering  phenomenon,  Ideally 
the  model  would  be  derived  directly  from  a  plasma  physics  theory  of  the  heater- 
ionosphere  interaction  and  would  be  generally  applicable  to  all  scattering 
situations  of  practical  interest.  It  will  probably  be  at  least  a  year  before 
the  relevant  theory  is  developed  to  the  point  where  this  can  be  done.  At  the 
otl  '  extreme,  the  model  could  be  given  simply  as  a  tabulation  of  experimental 
results.  The  problem  with  tiis  approach  is  that  data  is  available  over  only  a 
limited  number  of  paths  and  with  only  a  few  types  of  equipment.  A  tabulation 
of  experimental  results  provides  no  information  about  how  to  extrapolate  these 
results  to  other  situations.  As  an  interim  solution  we  have  divided  the 
modeling  task  into  three  parts, 

1 .  on-frequency  scattering 

2.  plasma  line  scattering 

3.  scintillation 

The  on-frequency  scattering  and  scintillation  will  be  handled  by  asking  what 
size,  shape  and  distribution  of  heater-produced  ionospheric  irregularities  will 
scatter  in  the  observed  way.  Given  these  irregularities,  the  on-frequency 
scatter  and  the  scintillation  under  other  experimental  conditions  can  be  pre¬ 
dicted.  The  theory  for  plasma  line  scattering  is  far  enough  along  so  that  in 
this  case  a  satisfactory  scattering  model  can  be  derived  directly  from  the 
theory.  This  paper  is  devoted  essentially  to  the  on-frequency  scattering  model 
while  the  models  for  the  latter  two  were  presented  by  Perkins1  and  Bowhill2, 
respectively. 
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II  PHYSICAL  DESCRIPTION  (U) 

(S)  The  model  has  been  developed  for  the  Platteville  transmitter  heating 
at  full  power  with  ordinary  ray  polarization  and  with  a  vertically  directed 
beam.  The  heater  frequency  is  assumed  to  be  at  or  below  the  critical  frequency 
of  the  ionosphere.  A  sketch  of  the  disturbed  volume  which  develops  under  these 
situations  is  shown  in  Figure  1.  For  the  sake  of  clarity,  the  heater  beam  has 
been  idealized  as  having  sharp  edges,  that  is,  the  gain  is  constant  within  the 
beam  but  drops  abruptly  to  zero  outside  of  the  beam.  The  heating  wave  re¬ 
flects  at  the  altitude  where  the  heater  frequency  equals  the  local  plasma  fre¬ 
quency.  The  irregularities  responsible  for  plasma  line  scatter  are 
generated  at  essentially  this  same  altitude  and  are  confined  there.  These 
irregularities  are  not  extremely  aspect  sensitive  and  consequently  to  a  radar 
the  plasma  line  scatter  appears  to  come  from  a  thin  disk  which  can  be  moved  up 
or  down  in  altitude  by  changing  the  heater  frequency.  A  second  group  of  irre¬ 
gularities  are  also  generated  at  the  heater  reflection  height  and  these  are 
responsible  for  on-*requency  scattering.  In  contrast  to  the  plasma-line  irre¬ 
gularities  (which  /.ire  probably  caused  by  propagating  waves)  the  on-frequency 
irregularities  are  essentially  stationary  in  space  and  extend  along  the  earth’s 
magnetic  field  for  several  tens  of  kilometers.  The  on-frequency  irregularities 
are  highly  aspect  sensitive  and  consequently  scatter  significantly  only  when 
specularity  is  achieved.  For  the  backscatter  case  shown  in  Figure  1,  this 
means  that  the  radar  must  look  perpendicular  to  the  magnetic  field  in  order  to 
detect  a  signal.  The  "ocus  of  points  for  which  the  radar  can  look  perpendi¬ 
cular  to  the  field  is  called  the  orthogonality  contour  and  this  is  fixed  in 
space  by  geography  alone,  having  nothing  to  do  with  the  heater  parameters.  The 
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(S)  FIGURE  1  SKETCH  OF  THE  RF  SCATTERING  MODEL  (U) 
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on-frequency  irregularities  are  most  intense  at  the  heater  reflection  height 
and  consequently  the  strongest  on-frequency  scatter  will  be  observed  when  the 
heater  reflection  level  is  adjusted  to  coincide  with  the  orthogonality  contour. 

Ill  ON-FREQUENCY  SCATTERING  MODEL  (U) 

(S)  The  on-frequency  scattering  model  has  been  arrived  at  by  applying 
the  Booker  theory3  for  scattering  from  anisotropic  irregularities  to  the  Ivory 
Coral  problem.  Booker  applied  his  theory  to  the  problem  of  radio  wave  scat¬ 
tering  from  the  aurora  and  found  that  the  experimental  data  could  be  under¬ 
stood  if  there  were  weak  irregularities  of  electron  density  present  that  were 
aligned  with  the  earth's  magnetic  field.  We  find  the  same  to  be  true  for  the 
heated  volume  over  Platteville.  As  with  the  auroral  case,  there  is  no  attempt 
made  to  explain  how  these  irregularities  are  produced  or  maintained.  This  is 

a  plasma  physics  question  and  is  being  addressed  by  the  Ivory  Coral  theorists 
at  the  present  time. 

(S)  The  scattering  model  we  present  is  a  volume  scattering  coefficient 
which  must  be  integrated  over  the  heated  volume  to  give  the  total  power  re¬ 
ceived  over  a  given  path.  Several  alternatives  exist  for  doing  the  required 
integration  in  practice.  For  special  cases  the  integration  can  be  done 
analytically  leading  to  simple  formulas  which  can  be  evaluated  on  a  slide  rule. 
For  example,  this  has  been  done  for  the  case  of  backscatter  measurements  made 
magnetically  south  of  the  heated  volume  (the  situation  under  which  most  of  our 
experimental  data  was  taken)  and  this  has  simplified  the  task  of  varying  the 
parameters  of  the  model  to  fit  the  backscatter  (or  near  backscatter)  observations. 
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It  is  likely  that  satisfactory  analytic  integrations  can  be  done  for  most 
specific  cases  of  interest  but  these  must  be  looked  at  one  at  a  time  to  decide 
what  idealizations  are  appropriate.  The  more  general  approach  is  to  use  a 
computer  to  numerically  integrate  the  volume  scattering  coefficient  over  the 
heated  volume.  In  this  way  realistic  antenna  patterns,  magnetic  field  models, 
and  ambient  electron  density  profiles  can  be  used.  If  this  approach  is 
followed  it  will  probably  be  important  to  minimize  the  computer  time  required. 
One  promising  way  of  doing  this  is  to  recognize  that  the  high  aspect  sensi¬ 
tivity  of  the  scatterers  means  that  only  that  fraction  of  the  heated  volume 
that  lies  near  the  specularity  surface  for  the  communications  link  being  studied 
need  be  considered  in  the  integral.  In  other  words,  the  volume  integral  can  be 
reduced  to  a  surface  integral  with  an  appreciable  saving  in  computer  time. 


(U)  The  volume  scattering  coefficient  (power  scattered  in  the  direction 
of  the  receiver,  per  unit  solid  angle,  per  unit  incident  power  density,  per 
unit  volume)  is^ 


(2tt) 


3/2  2 


(a  M/N) ‘ 


sin 


2 

a  ■pi 
1 1  L1 


8tt2s  i  n2 


exp- 


(T^cos^i  +L2s  i  n2i^ ) 


(1-«)T^L2 


8tt2s  i  n2-§- 


exp 


2  7  2  o  p  \ 

"  (T^cos  ip  +  L^sin  p) 


(1) 


where 

X  =  angle  between  the  direction  of  scattering  and  the  electric 
field  of  the  incident  wave 

e  =  angle  between  the  incident  and  scattered  wave  directions 
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(a  n/n) 
t,,t2 

L-|  ,  \-2 


a 


complement  of  the  angle  between  the  external  bisector  to 
incident  and  scattered  directions  and  the  magnetic  field 

operating  wavelength 

the  plasma  wavelength  corresponding  to  the  ambient  electron 
density 

the  mean  square  electron  density  fluctuations  of  the  irre¬ 
gularities 

Gaussian  autocorrel ation  scale  sizes  for  density  fluctuations 
transverse  to  the  magnetic  field 

Gaussian  autocorrelation  scale  sizes  for  density  fluctuations 
parallel  to  the  magnetic  field 

a  constant  determining  the  relative  strength  of  the  two 
Gaussian  irregularity  distributions  used  in  this  model 


(U)  There  are  three  angles  involved  in  the  model;  x.  0,  and  ip .  The  angle 
X  accounts  for  the  fact  that  the  irregularities  scatter  much  like  small  dipoles 
with  the  null  direction  along  the  E  field  of  the  incident  wave.  The  term  sin2x 
can  be  made  unity  by  suitable  choice  of  transmitter  and  receiver  polarizations. 
The  angles  e  and p  are  illustrated  in  Figure  2.  Let  the  scattering  plane  be 
defined  by  the  transmitter,  the  scattering  volume  element,  and  the  receiver. 
Since  this  is  weak  scattering,  most  of  the  incident  signal  continues  on  through 
the  scattering  volume  without  deviation  and  this  is  called  the  transmitted 
signal.  The  angle  between  the  transmitted  and  the  scattered  signal  is  defined 
as  e.  The  bisector  of  the  path  from  transmitter  to  scattering  volume  to  re¬ 
ceiver  also  lies  in  the  scattering  plane.  The  angle  between  this  bisector  and 
the  magnetic  field  is  defined  as  (90  )°. 
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'S)  FIGURE  2  BISTATIC  SCATTERING  GEOMETRY  (U) 
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(S)  The  operating  wavelength  of  the  communication  system  is  x.  The 
local  electron  density  is  given  in  terms  of  the  wavelength,  x^,  corresponding 
to  the  local  plasma  frequency.  XN  will  vary  through  the  heated  volume  and 
should  be  included  in  the  volume  integration. 


(S)  The  intensity  of  the  electron  density  fluctuations  is  given  by 


2 

(aN/N)  and  this  varies  throughout  the  heated  volume.  The  maximum  value, 

2 

(aN/N)  ,  occurs  at  the  heater  reflection  level  on  the  axis  of  the  heater  beam. 

At  points  off  the  axis  of  the  heater  beam  but  still  at  the  heater  reflection 

p 

level,  (a  N/N )  is  proportional  to  the  heater  power  density.  Thus  for  a  Gaussian 

p 

heater  beam,  (aN/N)  varies  in  a  Gaussian  manner  across  the  heater  reflection 
plane.  The  strength  of  t. <e  irregularities  drops  off  above  and  below  the  heater 
reflection  level  with  a  Gaussian  scale  height  A.  This  aspect  of  the  model  is 
illustrated  in  Figure  3. 


(S)  In  the  Booker  auroral  scattering  model  the  characteristics  of  the 
irregularities  are  specified  by  a  two-dimensional  autocorrelation  function 
having  a  correlation  length  T  transverse  to  the  earth's  magnetic  field  and  a 
con  elation  length  L  along  the  field.  Essentially,  T  controls  the  frequency 
dependence  and  L/T  controls  the  aspect  sensitivity  of  the  scatterers.  In  order 
to  match  the  frequency  dependence  observed  in  the  Ivory  Coral  phenomenon  we 
must  assume  that  two  sets  of  irregularities  are  present,  one  set  with  a  trans¬ 
verse  correlation  distance  T-j  and  another  with  a  transverse  correlation  dis¬ 
tance  T£.  The  experimental  measurements  to  date  give  no  evidence  for  aspect 
sensitivity  less  than  perfect.  That  is,  there  is  no  reason  to  assume  a  fi"-!t“ 

L  and  the  model  could  be  simplified  somewhat  by  using  a  one-dimens'  .  ?. jto- 
correlation  function  (variations  across  the  field  only).  However,  the 
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(S)  FIGURE  3  SPATIAL  VARIATION  OF  (AN/N)2  (U) 


SECRET 


41 


SECRET 


(S)  Cont'd. 

simplification  is  not  great  and  aspect  sensitivity  experiments  are  still  underway 
(HF-VHF  interferometry)  so  for  the  present  we  shall  retain  the  parameter  L  in 
the  model  and  assign  it  a  large  enough  numerical  value  to  produce  high  aspect 
sensitivity.  The  relative  strength  of  the  two  sets  of  Gaussian  irregularities 
in  this  model  is  set  by  the  parameter  a  . 

(S)  The  constants  (aN/N)^,  A,  T, ,  T2,  L, ,  L,  and  have  been  evaluated 
from  the  available  experimental  data  and  found  to  be  the  following: 

(aN/N)q  =  0.6  x  1  O'4 

A  =  10  kilometers 
T-|  =  0.1  meter 

T2  =  0.5  meters 
Lt  =  l_2  =  10  kilometers 
«  =  0.24 

The  bulk  of  the  quantitative  experimental  data  in  the  Ivory  Coral  program  has 
been  collected  over  backscatter  paths  or  over  bistatic  paths  which  are  so  close 
to  the  backscatter  geometry  that  the  difference  is  insignificant.  These  ex¬ 
perimental  results  are  summarized  by  the  shaded  band  (1)  shown  in  Figure  4.  The 
computed  backscatter  curve  using  the  on-frequency  scattering  model  with  the  con¬ 
stants  given  above  is  shown  by  the  solid  curve  labeled  (2).  The  model  predic¬ 
tion  is  seen  to  match  the  experimental  measurements  to  within  the  uncertainty 
of  the  data.  The  solid  curve  labeled  (3)  on  this  figure  is  the  model  prediction 
over  a  bistatic  path  for  which  the  scattering  angle  is  60°.  The  important  point 
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is  that  at  frequencies  such  as  200  MHz  the  radar  scattering  cross-section  over 
bi static  paths  is  expected  to  be  as  much  as  20  dB  greater  than  over  backscatter 
paths  at  the  same  frequency. 

(S)  The  cross-section  "nain"  expected  over  bistatic  paths  is  a  direct 
consequence  of  the  Bragg  scattering  law  and  the  observational  fact  that  the 
backscatter  cross-section  decreases  with  increasing  frequency.  The  Booker 
scattering  theory  views  the  irregularities  within  the  heated  volume  as  being 
a  super  position  of  plane  sinusoidal  waves  in  space,  only  one  component  of 
which  is  important  in  determining  how  strong  the  scatter  will  be  at  a  particular 
angle  and  wavelength.  The  important  component  is  that  which  satisfies  specularity 
for  the  path  (path  bisector  normal  to  scattering  surfaces)  and  which  has  a  spa¬ 
tial  wavelength  such  that  the  scattered  signal  from  each  cycle  adds  construc¬ 
tively  in  the  scattering  direction.  Constructive  interference  requires  that 
the  path  length  difference  for  waves  scattered  from  adjacent  planes  be  equal  to 
the  signal  wavelength.  Figure  5  shows  that  for  bistatic  paths  the  same  irre¬ 
gularity  component  that  supports  near-backscatter  for  an  operating  wavelength, 
xl»  wi11  support  scatter  over  a  more  bistatic  path  at  a  shorter  wavelength ,x2> 

The  upshot  of  this  is  that  the  observed  backscatter  curve  of  cross-section 
versus  frequency  can  be  used  to  predict  the  cross-section  curve  over  bistatic 
paths,  using  the  scaling  relation 

°e  ^  =  °180°  S1n  6/2) 

For  example,  if  e  =  60°,  sin  6/2  =  0.5  and  consequently  the  cross-section 
observed  at  frequency  f  over  a  bistatic  path  with  0  =  60°  will  be  the  same  as 
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that  observed  at  half  the  frequency  in  backscatter.  Sii.ce  the  cross-section 
increases  with  decreasing  frequency,  bistatic  paths  should  always  give  a  higher 
cross-section  than  backscatter  paths,  all  else  being  equal.  This  model  pre¬ 
diction  has  not  been  checked  experimental ly  but  plans  are  underway  to  do  so. 

IV  SUMMARY  AND  CONCLUSIONS  (U) 

(S)  An  on-frequency  scattering  model  has  been  developed  for  use  over 
bistatic  and  backscatter  paths.  The  model  is  consistent  with  the  experimental 
data  collected  over  backscatter  or  near-backscatter  paths.  Two  important  facets 
of  the  model  are  the  following: 

(S)  1.  The  aspect  sensitivity  is  great  enough  so  that  for  practical 

purposes  it  can  be  considered  perfect.  This  means,  for  example, 
that  the  width  of  the  bands  along  the  ground  where  signals  can 
be  received  from  a  remote  transmitter  via  the  heated  volume 
are  set  by  the  finite  dimensions  of  the  scattering  volume  rather 
than  by  the  finite  correlation  length  of  the  scatterers. 

(S)  2.  Significantly  larger  radar  cros. -sections  are  predicted  for 

bistatic  geometries  than  for  backscatter  geometries  at  the 
same  freo-jency.  For  example,  at  200  MHz  the  radar  cross-section 
of  the  heated  volume  should  increase  by  20  dB  going  from  back¬ 
scatter  (e  =  180°)  to  a  bistatic  path  (e  =  60°). 

(S)  The  model  is  felt  to  be  an  important  step  in  the  right  direction  but  is 

by  no  means  complete.  Except  for  special  cases,  a  numerical  integration  must  be 
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done  to  predict  signal  strengths  and  this  will  be  time  consuming  on  a  computer. 
This  drawback  of  the  model  can  probably  be  overcome  for  the  cases  of  interest 
by  simplifying  the  model  to  the  point  where  the  integration  can  be  done 
analytically.  More  fundamentally,  the  model  is  based  upon  the  assumed  exis¬ 
tence  of  irregularities  whose  characteristics  are  chosen  simply  to  match  the 
observations.  As  the  theory  of  the  heater-ionosphere  interaction  is  developed 
to  the  point  where  the  characteristics  of  the  irregularities  can  be  predicted, 
the  on-frequency  scattering  model  should  be  revised. 
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ABSTRACT  (U) 

(S)  The  paper  presents  a  method  to  extract  information  on  the  scale 
sizes  related  to  longitudinal  correlation  and  intensity  of  the  density  fluctua¬ 
tions  from  the  backscatter  radar  observations  of  the  angular  pattern  and  the 
orientation  of  the  reflection  surface.  The  two  parameters  determine  the  degree 
to  which  deviations  occur  in  the  angular  dimension  of  the  radar  pattern  and 
the  tilt  of  the  reflection  surface  from  that  expected  for  the  conditions  of 
perfect  aspect  sensitivity.  The  curves  generated  of  the  deviations  parametric 
in  the  two  scale  sizes  show  that  the  RAM  observations  are  consistent  with  a 
recently  reported  value  of  7.5  km  for  the  intensity  scale  size  and  suggest  a 
lower  limit  of  about  10  m  for  the  longitudinal  correlation  length  defined 
through  the  relation  6s  =  x/2  f^L. 
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I  INTRODUCTION  (U) 

(S)  The  'on-frequency'  scattering  model  presented  in  an  earlier  paper 
provided  estimates  for  the  transverse  scale  size  and  the  intensity  of  the 
density  fluctuations  for  the  heater  generated  irregularities.^  The  model  was 
based  on  the  assumption  that  the  signal  returns  are  extremely  aspect  sensitive 
and  hence,  it  turned  out  that  the  total  cross  section  is  independent  of  the 
longitudinal  scale  size.  The  intensity  of  the  scattering  irregularities  is 
believed  to  be  a  function  of  height  with  its  peak  located  at  the  height  of 
maximum  heating.  Although  this  dependence  is  built  into  the  mooel  ,  no  attempt 
has  been  made  to  obtain  a  scale  size  to  describe  the  intensity  of  the  density 
fl  uctuations . 

(S)  The  information  on  the  longitudinal  correlation  and  the  fluctuation 
intensity  scale  sizes  is  considered  to  be  essential  for  an  accurate  description 
of  the  scattering  model.  The  two  parameters  influence  in  an  important  manner 
the  intensity  and  the  angular  dimension  of  the  scattering  lobe  for  a  constant 
range  from  a  diagnostic  radar.  The  scatter  returns  will  originate  from  the 
region  where  the  radar  looks  transverse  to  the  field  aligned  irregularities  and 
the  scattering  lobe  will  have  the  dimension  of  the  antenna  pattern  for  the  con¬ 
ditions  of  perfect  aspect  sensitivity.  Any  deviation  from  this  exact  situation 
is  a  direct  consequence  of  the  finite  aspect  sensitivi  ly  which  in  turn  results 
from  the  finiteness  of  the  longitudinal  correlation  scale  size.  When  the  irre¬ 
gularities  are  highly  aspect  sensitive,  the  primary  effect  of  the  height  depen¬ 
dence  of  the  de.isity  fluctuations  is  to  cause  the  intensity  of  the  scattered 
signal  to  change  as  a  function  of  mismatch  between  the  heater  and  the  specular 
reflection  heights.  An  estimate  has  been  made  of  the  scale  size  for  the 
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fluctuation  intensity  by  Sweeney  and  Flaler  using  the  observations  of  the 
signal  intensity  variation  with  the  height  mismatch.2  A  somewhat  less  obvious 
effect  is  to  Introduce  a  slight  shift  In  the  angular  position  of  the  maximum 
return  in  the  direction  of  the  heater  reflection  height  from  the  point  of  exact 
specularity.  Thus  the  effects  of  the  two  parameters,  which  in  fact  are  coupled, 
will  reflect  in  the  deviations  in  the  orientation  of  the  reflection  surface  and 
the  angular  width  of  the  scattered  return  from  that  corresponds  to  the  condi¬ 
tion  of  perfect  aspect  sensitivity.  This  paper  is  devoted  to  the  simulation 
of  the  scattering  patterns  for  RAM  radar  and  to  extract  some  information  on 
the  two  parameters. 


II  SIMULATION  OF  RADAR  PATTERNS  (U) 

(S)  The  observed  angular  pattern  of  the  signal  intensity  is  the  result 
of  the  convolution  of  the  radar  beam  with  the  scattering  pattern  of  the  irre¬ 
gularities.  Hence  the  intensity  pattern  I(eo)  for  an  antenna  scan  parallel  to 
the  magnetic  field  (elevation  scan  for  RAM)  is  expressed  as: 

+  7T 

!(%)  =  K  j  G2  (e-8Q)  Gs  (9)  de  (1) 


where,  as  illustrated  in  Figure  1, 

Ga(e)  =  antenna  pattern  of  the  radar 

Gs(e)  =  scattering  pattern  of  the  irregularities 
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(S)  FIGURE  1.  ILLUSTRATION  TO  SHOW  THE  GEOMETRY  USED  IN 

PATTERN  AND  THE  SLOPE  OF  THE  SC 
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e0  =  elevation  angle  of  the  radar  beam  from  the  direction 
of  specularity 

e  =  elevation  angle  to  a  scattering  irregularity 
and  K  ~  a  constant  related  t.  system  parameters 

Let  the  radar  beam  and  the  scattering  patterns  be  of  Gaussian  shape  such  that : 

Ga  =  <jo  exp  {-<6-9o>2'ea  }  (2) 

and 

G$  =  G$(y)  exp  j  -(2e/es)2  }  (3) 

Where  efl  and  es  are  the  e-folding  beamwidths  of  the  antenna  and  the  scattering 
patterns.  G$(y)  is  the  scattering  lobe  intensity  corresponding  to  the  intensity 
variation  of  the  irregularities  as  function  distance  along  the  magnetic  field. 
Let  us  assume 

Gs <y>  =  Gso  exp  { -(y-y0)2/*2}  (4) 

where,  as  shown  in  Figure  1, 

yQ  =  distance  of  the  heater  reflection  from  the  specular  point 

y  =  distance  from  the  specular  point  to  the  scattering  irre¬ 
gularity 

and  A  =  Gaussian  scale  size  for  the  intensity  variation  of  the 

irregularities  along  the  magnetic  field 
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The  high  aspect  sens  tivlty  of  the  irregularities,  limiting  the  scattering  to 
a  spatial  extent  small  compared  to  the  radar  range,  permits  the  approximation 
e  *y/R  for  evaluating  the  integral  in  Equation  (1).  On  substituting  the  above 

relations  and  evaluating  the  integral  one  obtains  the  following  for  the  signal 
intensity  pattern: 


pRe 

j— 

LS(R9a)‘ 


(5) 


where 


1/2 


The  angular  position  6m  of  the  peak  of  the  intensity  pattern  is  obtained  by 

solving  for  9o  satisfying  the  equation  dl'e  )/de  =  0.  The  solution  for  e  is 

.  U  0  m 

given  as: 


em  =  2%ea2/A2[2S  2(Rea)2  -  4]  (6) 

The  intensity  distribution  is  expressed  in  the  normalized  form  as: 

Ve0>  =  He^/Uem)  (7) 

using  Equations  (5)  and  (7)  the  normalized  intensity  patterns  can  be  computed 
parametric  in  A  and  the  scale  sizes  for  the  fluctuation  intensity  and  the 
longitudinal  correlation  of  the  irregularities.  The  correlation  length  L  is 
related  to  the  scattering  lobe  size  eg  as^: 

es  =  */2V2VL  (8) 
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Where  x  is  the  wavelength  of  the  operating  radar.  Figure  2  presents  a  set  of 
curves  parametric  in  L  showing  the  variation  of  es  with  the  operating  frequency. 
The  angular  dimension  and  the  slope  of  the  scattering  disk  of  Irregularities 
are  determined  from  the  computed  normalized  intensity  patterns  I  (e  ).  The 

of  the  scattering  disk  is  determined  from  the  angular  position  em  of  the 
peak  of  the  intensty  pattern. 

Ill  NUMERICAL  RESULTS  (U) 

(S)  The  following  parameters  relevant  to  the  RAM  radar  observations  are 
used  to  compute  the  intensity  patterns: 

R  =  900  km 
ea  =  3.46° 
f  =  157.5  MHz 
yQ  =  -10  km 

The  patterns  are  calculated  parametric  in  the  scale  sizes  A  and  L.  The  computed 
values  of  the  angular  width  and  the  tilt  angle  of  the  scattering  disk  are 
normalized  by  the  antenna  beamwidth  and  the  tilt  of  the  specular  scattering 
surface  and  the  results  are  presented  in  Figure  3.  The  curves  generated 
correspond  to  the  scale  sizes  A  ranging  from  5  to  100  km  and  L  from  2.5  to  16  m. 
The  dependence  of  the  width  and  the  tilt  of  the  scattering  disk  on  the  scale 
size  parameters  is  found  to  be  sensitive  enough  to  extract  information  on  the 
scale  sizes  from  the  backscatter  radar  measurements.  It  is  recognized,  however, 
that  at  very  high  operating  frequencies  the  observed  deviations  from  the  condi¬ 
tions  of  perfect  aspect  sensitivity  are  very  small  and  the  method  relies  on  a 
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(U)  FIGURE  2.  FREQUENCY  DEPENDENCY  OF  THE  ASPECT  SENSITIVITY  OF  THE 

SCATTERING  IRREGULARITIES  PARAMETRIC  TO  THE 
LONGITUDINAL  AUTOCORRELATION  LENGTH  L  (U) 
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N-S  SLOPE  OF  SCATTERING  DISK 
N-S  SLOPE  OF  ORTHOGONALITY  SURFACE 


(s)  FIGURE  3.  THE  NORMALIZED  ANGULAR  DIMENSION  AND  THE  SLOPE 
OF  THE  SCATTERING  DISK  PARAMETRIC  IN  THE  SCALE 
SIZES  A  AND  L  (S) 
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high  degree  of  precision  in  the  measurements.  The  received  power  pattern  in 
elevation  scan  observed  by  Minkoff  et.  al4  using  the  RAM  radar  during  a  Prairie 
Smoke  II  test  conducted  on  April  25,  1972  is  shown  in  Figure  4.  The  circles 
represent  the  received  power  and  the  crosses  denote  the  square  of  the  antenna 
pattern  corresponding  to  a  point  source  response  for  the  backscatter  operation. 
The  dotted  curves  represent  10%  and  20%  broadening  over  the  antenna  pattern. 

The  observations  do  suggest  that  the  broadening  due  to  finite  aspect  sensitivity 
of  the  scatterers  is  almost  negligible  and  may  be  not  more  than  5%  at  the  most. 
The  results  on  the  tilt  angle  of  the  scattering  disk  are  taken  from  the  obser¬ 
vations  collected  at  RAM  on  12,  13  and  14  of  October  1971  . 5  The  observed  tilt 

angles  normalized  to  the  tilt  angle  of  the  specular  surface  are  shown  in  Figure 

5  for  nine  different  cases.  It  may  be  seen  from  the  figure  that  the  deviation 

of  the  observed  tilt  angle  relative  to  the  specular  surface  is  15%  or  less. 

The  observations  are  found  to  be  consistent  with  a  recently  reported  estimate 
of  7.5  km  for  A  in  that  it  requires  only  a  little  over  1%  broadening  in  the 
measured  pattern.  Commensurate  with  these  parameters  and  a  15%  deviation  in  the 

tilt  angle,  the  curves  suggest  a  lower  limit  of  about  10  m  for  the  correlation 
scale  size  L. 


IV  CONCLUSION  (U) 

(S)  It  is  shown  that  the  backscatter  radar  observations  of  the  angular 
dimension  of  the  radar  pattern  and  the  tilt  of  the  reflection  surface  can  be 
used  to  extract  information  on  the  characteristic  scale  sizes  describing  the 
longitudinal  correlation  and  the  intensity  variation  of  the  density  fluctuations 
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FREQUENCY  OF  OCCURRENCE 


FIGURE  4.  RECEIVED  POWER  VERSUS  ELEVATION  ANGLE  OBSERVED  AT  RAM  ON  25  APRIL 
1972  ALONG  WITH  CURVES  SHOWING  10*  and  20%  PATTERN  BROADENING  (S) 


N-S  SLOPE  OF  ORTHOGONALITY  SURFACE 


(Sy  FIGURE  5.  DISTRIBUTION  OF  THE  NORMALIZED  SLOPE  OF  THE  SCATTER  INC  DTSk: 

OBSERVED  AT  RAM  ON  12,  13  AND  14  OCTOBER  19T  (S) 
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RAM  observations  are  found  to  be  consistent  with  a  value  of  7.5  km  for  the 
intensity  scale  size  and  seem  to  suggest  a  lower  limit  of  10  m  for  the  longi¬ 
tudinal  correlation  length  as  defined  by  Booker.3 
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DIMENSIONS  OF  GROUND  ILLUMINATION  PATTERNS  (U) 

by 

W.  H.  Russell,  B.  R.  Pendyala  and  G.  D.  Thome 

Raytheon  Company 
Sudbury,  Massachusetts 

ABSTRACT 

(S)  This  paper  presents  a  method  for  computing  the  ground  illumination 
patterns  for  'on-frequency'  scatter  from  the  heated  volume.  The  size  of  the 
illumination  zones  in  a  north-south  plane  are  calculated  by  considering  their 
dependency  on  the  operating  frequency,  the  size  of  the  scattering  volume  and 
the  aspect  sensitivity  of  the  scatterers.  The  results  show  that  the  size  of 
the  illumination  zones  are  primarily  dependent  on  the  finite  size  of  the 
scattering  volume  rather  than  on  the  finite  aspect  sensitivity  of  the  scatterers. 
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I  INTRODUCTION  (U) 

(S)  Aircraft  observations  of  scatter  from  the  heated  region  during  June 
of  1971  have  shown  that  adequate  reception  of  the  scattered  energy  is  possible 
only  within  specific  zones  on  the  ground.^  These  regions  are  typically  200  to 
250  km  wide  in  the  north-south  direction  but  considerably  larger  in  the  east- 
west  direction.  It  was  postulated  at  that  time  that  the  size  of  the  north- 
south  zone  was  dependent  primarily  on  the  size  of  the  scattering  volume.2 
Experimental  evidence  also  showed  that  operating  frequency  influenced  the  zone 
size  and  that  refraction  apparently  affects  the  location  of  the  zones.  Other 
factors  such  as  the  aspect  sensitivity  and  the  propagation  geometry  will  also 
influence  the  size  of  the  illumination  area  and  its  location. 

II  COMPUTATION  OF  GROUND  ILLUMINATION  PATTERNS  (U) 

(S)  To  investigate  some  of  these  interdependencies,  the  scattering  pro¬ 
cess  has  been  modeled  analytically  so  that  the  relative  intensity  of  the 
scatter  on  the  ground  could  be  evaluated.  To  examine  only  the  size  of  the 
zone  and  how  it  varies  with  operating  frequency  and  the  size  of  the  heated 
volume,  some  approximations  have  been  made.  The  frequency  dependency  has  been 
restricted  to  the  VHF  spectrum  so  that  refraction  effects  can  be  neglected. 

Also  the  scattering  geometry  has  been  constrained  to  a  plane  containing  the 
heated  region  and  the  radar.  These  approximations  permit  a  two  dimensional 
straight  line  description  of  the  scattering  process.  Furthermore,  the  magnetic 
field  has  been  considered  constant  throughout  the  heated  volume  and  the  curva¬ 
ture  of  the  earth  was  neglected. 
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(s)  Within  the  above  constraints,  the  resulting  geometry  of  the  scat¬ 
tering  may  be  described  by  Figure  1.  The  intensity  of  heater  produced  irre¬ 
gularities  is  considered  maximum  at  point  0  and  decreases  in  a  Gaussian  fashion 
in  the  x  and  z  direction  with  scale  sizes  of  W  and  A  respectively. 

(S)  The  intensity,  A  I,  at  a  point  Q  on  the  ground  due  to  a  scatterer  at 
P  is,  therefore,  dependent  on  the  location  of  P  in  the  scattering  volume  so 
that  A  I  may  be  set  proportional  to  the  following  term,  (see  Figure  1) 

-(x2/W2)  -  (z  +  x  cos  I)2/ A2 

6  (1) 

The  intensity  will  also  be  influenced  by  the  scattering  pattern  of  the  irre¬ 
gularities  at  P.  For  a  Gaussian  scattering  pattern  of  the  irregularities  with 
a  1/e  beamwidth  of  e$ ,  the  effect  of  aspect  sensitivity  on  the  relative  inten¬ 
sity  at  Q  may  be  described  by  an  expression  of  the  form 


R  ~  R0  +  (x  -  d)  sin  I  (4) 

The  relations  in  Equations  3  and  4  are  appropriate  for  the  case  of  n^'gh  aspect 
sensitivity  and  when  it  is  assumed  that  Rq  is  long  compared  to  the  scattering 
volume  size  and  the  ground  distance  d. 
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(S)  From  the  expressions  1  and  2  the  intensity  at  Q  due  to  the  scatterer 


at  P  is  given  as: 


A  I  =  I  e 
m 


-  x2/W2  -  (z  +  x  cos  I )2/A2  -  (S/Re  )2 


Evaluating  the  above  expression  for  various  values  of  d  provides  a  distribution 
of  the  intensity  on  the  ground  due  to  scatter  from  a  single  irregularity  at  P. 
The  resulting  distribution  will  not  only  be  dependent  on  the  geometry  but  also 
on  the  operating  frequency  through  the  aspect  sensitivity  term  e  since3: 

e$  =  a/2\/?ttL  (6) 

where 

A  =  operating  wavelength 

L  =  autocorrelation  length  as  defined  by  Booker 

A  distribution  in  d  may  also  be  obtained  for  the  scatter  due  to  a  column  of 
irregularities  along  the  field  lines  by  integrating  Equation  5  over  all  z  so 
that 

ffIm  e  N  -x2/W2  -x2cos2I/A2  -d2cos2I/e  R2 


where 


m  d  sin  I 

M  ■  “2  +  2  -  -  /e?R2 

/  A  Rq  +  X  sin  I  s 


x  cos  I 

N  =  - —  + 

fc 


Yz---—1  ) 

A  Rq  +  x  sin  1/ 


d  cos  I  /e2R2  /M 
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(S)  In  a  similar  manner,  a  distribution  of  the  intensity  on  the  ground 
due  to  scatter  from  a  heated  area  may  be  obtained  by  integrating  the  effect  of 
scatterer  along  z  and  also  x.  The  integration  in  the  x  dimension,  however,  has 
been  accomplished  by  a  numerical  technique. 

(S)  To  illustrate  this  approach,  the  intensity  distributions  have  been 
computed  for  a  geometry  appropriate  to  that  of  the  RAM backscatter  observations 
of  the  irregularities  over  the  Platteville  heater.  Computations  of  the  dis¬ 
tributions  were  made  at  50,  100  and  150  MHz  for  a  single  irregularity  at  the 
center  of  the  heated  region,  a  column  of  irregularities  through  the  center  of 
the  heated  volume,  and  for  a  total  area  of  scatterers  centered  over  the  heater. 
The  parameters  chosen  to  compute  these  distributions  are  as  follows: 

A  =  7.5  km  R  =  900  km 

W  =  50  km  i  =  68° 

The  results  are  shown  in  Figure  2,  where  the  areas  under  the  distribution  curves 
are  shaded  to  distinguish  between  the  three  conditions  considered  above  (solid  = 
single  irregularity,  dark  dots  =  column,  and  light  dots  *  total  area).  It  is 
apparent  from  Figure  2  that  consideration  of  a  large  scattering  area  produces 
a  substantial  distribution  on  the  ground.  The  distributions  generated  from  the 
integration  over  the  total  heated  area  are  approximately  235  km  wide  between 
b  dB  points  which  compares  favorably  with  the  experimental  observations  of  zone 
sizes  of  approximately  200  to  250  km.  It  may  also  be  noted  from  Figure  2  that 
these  distributions  are  all  approximately  of  the  same  size  regardless  of  fre¬ 
quency.  This  condition  exists  since  any  broadening  of  the  distributions  due 
to  finite  highly  aspect  sensitive  scattering  patterns  of  the  individual 
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irregularities  is  small  compared  with  the  effect  of  the  total  scattering  volume. 

Ill  CONCLUSIONS  (U) 

(S)  The  results  presented  here  suggest  that  the  size  of  the  ground 
illumination  pattern  for  on-frequency  scattering  from  the  heated  volume  is  con¬ 
trolled  primarily  by  the  finite  size  of  the  scattering  volume  rather  than  by 
the  finite  aspect  sensitivity  of  the  individual  scatterers. 
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RF  SCATTERING  MODEL  SUMMARY  AND  CONCLUSIONS  (U) 

by 
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Sudbury,  Massachusetts 

(S)  An  RF  scattering  model  which  is  consistent  with  the  existing  experi¬ 
mental  measurements  has  been  developed  for  the  'on-frequency'  signal  component 
scattered  from  the  heater  induced  ionospheric  irregularities.  The  model  is 
intended  primarily  for  scatter  communication  systems  application  and  further 
to  aid  in  the  theoretical  and  experimental  efforts  to  understand  the  heater- 
ionosphere  interaction  phenomenon. 

(S)  The  model  presented  here  is  based  on  the  theory  of  weak  scattering 
(Born  approximation)  developed  by  Booker  for  the  case  of  anisotropic  field 
aligned  Gaussian  irregularities.  The  same  type  of  functional  form  has  been 
adopted  in  the  model  to  characterize  the  heater  induced  'on-frequency 'scat¬ 
tering  irregularities.  No  attempt  has  been  made  to  justify  how  such  type  of 
irregularities  are  produced  or  maintained  and  to  this  extent  it  might  become 
necessary  to  modify  the  model  in  future  when  a  complete  theoretical  descrip¬ 
tion  of  the  irregularities  is  made  available. 

(S)  The  characteristics  of  the  irregularities  in  the  heated  volume  are 
specified  generally  by  a  two  dimensional  Gaussian  autocorrelation  function  having 
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a  scale  size  T  transverse  to  the  magnetic  field  and  a  scale  size  L  along  the 
field.  These  two  parameters  along  with  the  intensity  of  the  density  fluctua¬ 
tions  determine  the  scattering  properties  of  the  disturbance.  Essentially,  it 
is  the  transverse  scale  size  that  governs  the  frequency  dependence  whereas 
the  ratio  (L/T)  determines  the  aspect  sensitivity  of  the  scatterers.  In  order 
to  be  able  to  matci  the  model  computations  with  the  observations  of  total 
scattering  cross  section  as  function  of  radar  frequency,  it  was  round  essential 
to  assume  that  two  dominant  sets  of  irregularities  with  transverse  scale  sizes 
T1  and  T2  are  Present  in  the  scattP“ing  volume.  Using  a  highly  refined  earth's 
magnetic  field  model  and  the  observations  of  the  distribution  of  signal  returns 
r  space,  it  was  shown  that  the  on-frequency  scatterers  are  extremely  aspect 
sensitive.  This  property  has  been  used  to  effect  substantial  simplification 
in  analytical  evaluation  of  the  scattering  coefficient  integral.  The  total 
backscatter  radar  cross  soction  for  two  scale  size  scattering  medium  was 
derived  as: 


4tt5W2IU 


(an/nV 


o  T  = 


exp 


{-dV} 

o-«) 


+ 


; 


where 

R  =  range  from  radar  to  the  specular  point  in  the  scattering 
volume 

a  -  radar  wavelength 
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W 


(aN/N)* 


A 

D 


a 


ambient  plasma  wavelength  at  the  specular  point 

Gaussian  r?dius  of  the  heated  volume  in  a  horizontal  plane 
at  the  reflection  level 

mean  square  electron  density  fluctuations  of  the  irregularities 
at  the  specular  point 

Gaussian  scale  size  of  (aN/N)^  along  the  magnetic  field 

separation  along  the  field  between  the  specular  point  and 
the  heater  reflection  level 

Gaussian  autocorrelation  scale  sizes  for  density  fluctuations 
transverse  to  the  magnetic  field 

a  parameter  determining  the  relative  strength  of  the  two 
Gaussian  irregularity  distributions  used  in  the  model 


(S)  The  model  predictions  were  found  to  match  well  with  the  observed 
optimum  radar  cross  sections  when  the  values  assigned  to  various  parameters  are 
as  follows: 


R  =  900  km 
=  54  m 
W  =  90  km 
D  =  0 


(AN/N)^  =  0.6  x  10"4 
T-j  =  0.1  m 

T 2  -  0.5  m 
a  =  0.24 


(S)  The  model  discussion  as  indicated  above  was  limited  essentially  to 
the  backscatter  case  although  an  approach  to  obtain  a  general  solution  for  a 
bistatic  situation  was  also  presented.  While  the  exact  solution  is  quite 
involved  it  can  be  shown  using  the  principle  of  Bragg  scattering  that  to  the 
first  order  of  approximation  the  radar  cross  section  for  a  bistatic  path  can 
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be  expressed  in  terms  of  an  equivalent  backscatter  cross  section  using  the 
scaling  relation 


CTe  Rp  R2)  =  a18Q  (f  sin  0/2,  R) 

Where  e  is  the  scattering  angle  and  R  is  the  equivalent  backscatter  range  related 
to  the  bistatic  ranges  R]  and  R£  as: 

R  ■  2  R]  R2  sin  (e/2)/(R]  +  Rg) 

The  scaling  relation  predicts  that  the  bistatic  radar  cross  sections  are  signi¬ 
ficantly  greater  than  that  for  backscatter  at  the  same  frequency  as  a  consequence 
of  the  fact  that  the  radar  cross  section  increases  with  decreasing  frequency. 

For  example,  at  200  MHz  the  radar  cross  section  of  the  heated  volume  should 
increase  by  approximately  20  db  going  from  backscatter  (e  =  180°)  to  a  bistatic 

path  (e  =  60°).  This  model  prediction  is  now  in  the  process  of  being  checked 
experimental ly. 


(S)  The  model  as  described  above  has  not  provided  any  information  on  two 
of  the  scale  size  parameters,  L  and  A,  describing  the  longitudinal  correlation 
and  the  intensity  variation  of  the  density  fluctuations.  A  method  was  pre¬ 
sented  to  extract  the  scale  sizes  from  the  backscatter  radar  observations  of 
the  angular  pattern  and  the  orientation  of  the  scattering  surface.  The  two 
parameters  determine  the  degree  to  which  deviations  occur  in  the  observed  radar 
patterns  from  that  expected  for  the  conditions  of  perfect  aspect  sensitivity. 
The  curves  generated  of  the  deviations  parametric  in  the  two  scale  sizes  show 
that  the  RAM  radar  observations  are  consistent  with  a  recently  reported  value 
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of  7.5  km  for  the  intensity  scale  size  and  suggest  a  lower  limit  10  m  for  the 
longitudinal  correlation  length  defined  through  the  relation  es  =  x  /2f?VL. 

(S)  The  final  aspect  of  the  scattering  model  presented  in  this  report  is 
related  to  the  dimensions  of  the  ground  illumination  patterns  due  to  the  'on- 
frequency'  scatter  from  the  heated  volume.  The  sizes  of  the  illumination  zones 
in  a  north-south  plane  are  calculated  by  considering  their  dependency  on  the 
operating  frequency,  the  size  of  the  scattering  volume  and  the  aspect  sensi¬ 
tivity  of  the  scatterers.  The  results  show  that  the  size  of  the  illumination 
zones  are  primarily  dependent  on  the  finite  size  of  the  scattering  volume 
rather  than  on  the  finite  aspect  sensitivity  of  the  scatterers. 
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I  1  3  AB5TR  A C  T 


An  scattering  model  which  is  consistent  with 
IVORY  CORAL  and  PRAIRIE  SMOKE  experimental  measurements  has  been 
developed  for  the  'on  frequency'  signal  component  scattered  from  the 
heater  induced  ionospheric  irregularities.  The  model  is  intended  pri- 
marily  for  scatter  communications  systems  application  and  further  to 
id  m  the  theoretical  and  experimental  efforts  to  understand  the 
heater- ionsphere  interaction  phenomenon. 

(S)  The  model  presented  here  is  bcsed  on  the  theory  of  weak 
scattering  (Born  approximation)  developed  by  Booker  for  the  case  of 
anisotropic  field  aligned  Goussian  irregularities.  The  characteristics 
of  the  irregularities  in  the  heated  volume  are  specified  generally  by 
a  two  dimensional  Gaussian  autocorrelation  function  having  a  scale  size 
transverse  to  the  magnetic  field  and  a  scale  size  L  along  the  field. 

ihe  model  predictions  were  found  to  match  well  with  the  observed  radar 
cross-sections. 


SECRET 

Security  Classification 


